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Zusammenfassung 



Die Anwesenheit verschrankter Quantcnzustande fiihrte zu Unbehagen und 
Zwcifcln an der Quantcntheorie, ob der Nichtlokalitat die sic in sich tragen. 
Da Lokalitat einer der Hauptpfciler physikalischer Gedankcngebaudc war und 
noch ist, wurde die Annahmc der UnvoUstandigkeit der Quantenthcoric und 
damit die Existenz versteckter lokaler Variablen geboren. Es war die Bahnbrc^- 
chende Arbeit von John Bell, die es moglich machte, anliand von sogonann- 
ten „Bellungleichungen" die Vorhersagen der Quantentheorie von denen einer 
Theoric mit versteckten Variablen zu unterscheiden; bislang sind derartige Mcs- 
sungcn zugunstcn dor Quantentheorie ausgefallen. Es ist seit der Griindung 
der Quanteninformationstheorie, dafi die Verschranktheit von Zustanden - das 
„ Entanglement" - wieder in den Mittelpunkt des Interesses geriickt ist, und 
zwar als Ressource fiir die Ausfiihrung klassisch unmoglicher Prozesse, wie et- 
wa der Teleportation. Dieser Aspekt fiihrte unausweichlich zu dem Bestreben 
wissenschaftlicher Arbeiten, diese Ressource klassifizicren und natiirlich auch 
quantifizieren zu woUcn. Zu diesem Zwecke sind Minimalkritericn an cin Mafi 
fiir Entanglement crarbeitet worden, welche auf das Konzept des „ Entanglement 
Monotones" fiihrten. Dieses Fundament motivierte nachfolgend lebhafte wissen- 
schaftliche Aktivitat mit Fokus auf das Entanglement von Bipartitionen, welche 
wichtige Kriterien wie den Schmidt-Rang, die von Neumann Entropie und die 
„ Concurrence" hervorbrachten. Der Erfolg bei bipartiten Systemen rief nach 
Erweiterungen der dort gefundenen Resultate auf multipartite Systeme; jedoch 
erwies sich dieses neue Feld als weitaus komplizierter. Ausschlaggebend hierfiir 
ist letztlich die parallele Existenz verschiedener Entangelmentklassen beziiglich 
stochastisch lokaler Transformationen begleitet von klassischer Kommunikation 
(SLOCC). 

Von vor etwa zehn Jahren entstammte die Idee, dal3 die Quanteninformati- 
onstheorie das Potential haben konnte, ein tieferes Verstandnis von komplexen 
Phanomenen im Bereich der kondensierten Materie oder der Quantenfeldtheorie 
zu erlangen. Tatsachlich fiihrte die darauffolgende intensive Forschungsarbeit 
unter dieser Pramisse auf beiden Gebieten zu einer wechselseitigen Befruch- 
tung. Von der dramatisch anwachsenden Intensitat wissenschaftlicher Arbeit auf 
dem Uberlapp beider Gebiete profitierten beide Seiten. Als besonders relevant 
fiir die vorliegende Arbeit sei hier die Untersuchung von Entanglementaspek- 
ten in der Nahe quantenkritischer Punkte zu nennen; jedoch fiihrte die Sicht- 
weise der Quanteninformationstheorie auch schon zu wichtigen Modifikationen 
numerischer Simulationsmethoden im Bereiche der kondensierten Materie, wie 
beispielsweise der DMRG. Desweiteren ist der vielerseits ertraumte Quanten- 
rechner letztendlich ein grol3es System von Quanteninformationseinheiten (z.B. 
Qubits), fiir welche lokale Operationen, aber auch paarweise Wechselwirkungen 
untereinander auf kontroUierte Weisc manipuliert worden konnen miissen. Daher 
sind z.B. Spinketten als Quantenregister, also als Trager von Quanteninforma- 
tion, bzw. als Quantenkanal vorgeschlagen worden. Im letzteren Falle wiirde die 
natiirlich gegebene hamilton'sche Zeitentwicklung zum Transport von Quanten- 
bits ausgenutzt werden woUen. 

Viele Arbeiten untersuchten also die Dynamik von Entanglement in Syste- 
men kondensierter Materie z.B. unter dem Aspekt optimaler Dateniibertragung, 
oder aber des Entanglementgehalts von Grundzustanden quantenkritischer Mo- 
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delle. Ein wichtiges und praktisch einhelliges Resultat der Untersuchungen der 

Ictztercn Katcgoric ist, dafi das fiir Quantcnphascniibcrgangc wichtigc Entan- 
glement vornehmlich multipartiter Natur ist. Diese Erkenntnis entfesselt eine 
Plethora ungeloster Probleme, welche bis in Bereiche der Invariantentheorie rei- 
chcn. Olinc klarc Vorstcllung, wclchc Entanglemcntklasse fiir bestimmte komple- 
xe Phanomene von Wichtigkeit sein konnte, wird ein moglicher Zusammenhang 
nur schwer hergestellt werden konnen; aber dafiir ware eine bekannte Klas- 
sifizierung dcs Entanglements Voraussetzung. Dieses Manko fiihrt dazu, dafi 
zunachst leicht berechenbare Grofien analysiert werden, welche aber dennoch 
gewisse Schliisse iiber das Entanglement im betrachteten System zulassen. Die- 
se erzwungen pragmatische Herangehensweise ist zwar wichtig; sie lafit jcdoch 
viele Facetten des Problems aus, und diirfte daher auf lange Sicht unzureichend 
sein. 

Die vorliegende Arbeit ist eine Zusammcnfassimg der in der vorhcrgehen- 
den Erklarung enthaltenen Liste von Publikationen. Nach einer Einfiihrung in 
die meistgenutzten Entanglementmafie, greift sie an das Problem des mutmafi- 
lichen Zusammenhangs zwischen Entanglement und Quantenphaseniibergangen 
an und zitiert eine Reihe von Arbeiten zu diesem Thema als Beleg fiir die Rele- 
vanz multipartiten Entanglements. Darauffolgend wird das Problem der Quan- 
tifizierung und Klassifizierung ,,genuin multipartiten Entanglements" formuliert 
und angegangen. Die Schliisselerkenntnis hierzu ist, die SL{2, C) samt Qubitper- 
mutationen als Invarianzgruppe zu identifizieren. Lokale antilineare Operatoren 
mit verschwindenden Erwartungswerten auf dem gesamten lokalen Hilbertraum 
werden als Bausteine fiir solche Mafie vorgestellt. Auf diese Weise konnte ei- 
ne voUstandige Klassifizierung vierpartiten Entanglements erfolgen, aber auch 
Mafic fiir echt multipartitcs Entanglement fiir cine bclicbigc Anzahl von Qubits 
sind in Reichweitc. Die konstruicrtcn Mafic sind zunachst nur wohldcfiniert auf 
reinen Zustanden; die Erweiterung auf gemischte Zustandc mittcls dcs sogenann- 
ten „convex roof" stellt eines der ungclostcn Probleme dar. Auf dem Wege zu 
dessen Losung konnten Gemische zweier bestimmtcr tripartiter Zustandc analy- 
tisch behandelt werden; einige der Ergcbnissc dicscr Arbeit lassen sich sogar auf 
belicbigc Rang-zwci Gemische und fiir bclicbigc Anzahl der Qubits iibcrtragen. 

Die Konstruktion mit lokalen invarianten Operatoren ist vom physikalischen 
Standpunkt bcsondcrs sinnvoU, da sie den Grunclstcin dafiir legt, die Entangle- 
mentmafie durch Korrclationsfunktionen auszudriickcn. Dazu konnte cine cins- 
zu-eins-Bezichung von Erwartungswerten antilinearer hcrmitcschcr Operatoren 
mit Erwartungswerten eindcutig zugcordneter linearer hcrmitcschcr Operato- 
ren hergestellt werden. Abschatzungen fiir das convex roof schwach gemischter 
Zustandc werden dann eine direkte Ankniipfung an das Experiment liefern. 
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Kapitel 1 

Introduction 



Entanglement has been always considered as the counter intuitive part and 
"spooky" non- locality inherent to quantum mechanics |15j . It apparently con- 
tradicts one of the basic pillars of physics - locality - and gave rise to severe 
skepticism for several decades. It was after the seminal contribution of John 
Bell that the fundamental questions related to the existence of entangled states 
could be tested experimentally. Under fairly general assumptions, Bell derived a 
set of inequalities for correlated measurements of two physical observables that 
any local theory should obey. The overwhelming majority of experiments do- 
ne so far are in favor of quantum mechanics thus demonstrating that quantum 
entanglement is physical reality [97^] . 

Entanglement has gained renewed interest with the development of quantum 
information science f5B]. In its framework, quantum entanglement is viewed at 
as a precious resource for quantum information processing. It is e.g. believed to 
be the main ingredient to the quantum speed-up in quantum computation and 
communication. Moreover several quantum protocols, as teleportation |18j . can 
be realized exclusively with the help of entangled states. 

The role of entanglement as a resource in quantum information has stimula- 
ted intensive research that tries to unveil both its qualitative and quantitative 
aspects [Ml [111 im I146| I100| llOlj . Many criteria have been proposed to distin- 
guish whether a pure state is entangled or not, as for example the Schmidt rank 
and the von Neumann entropy, and necessary requirements to be satisfied by an 
entanglement measure have been elaborated and have lead to the notion of an 
entanglement monotone |133| . Since then, a substantial bulk of work appeared 
on entanglement monotones for bipartite systems, in particular for the case of 
qubits. The success in the bipartite case for qubits asked for extensions to the 
multipartite case, but the situation proved to be far more complicated: different 
classes of entanglement occur, which are inequivalent not only under determini- 
stic local operations and classical communication (LOCC), but even under their 
stochastic analogue (SLOCC) [20] . 

During the last decade it has been suggested that quantum information 
science might bear the potential to give further insight into areas of physics as 
condensed matter or quantum field theory |104j . Indeed has the growing inte- 
rest of the quantum information community in systems from condensed matter 

^ There are states that do not violate Bell inequalities and nevertheless are entangled | 81| . 
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stimulated an exciting cross-fertilization between the two areas; the amount 
of work at the interface between condensed matter physics and quantum in- 
formation theory has grown tremendously during the last few years, shining 
light on many different aspects of both subjects. Methods developed in quan- 
tum information have proved to be useful also in the analysis of many-body 
systems. At the same time, the experience built up over the years in condensed 
matter physics is valuable for finding new protocols for quantum computation 
and communication: at the end, a quantum computer will be a many-body sy- 
stem where, differently from the 'traditional' situation, the Hamiltonian must 
be manipulated in a controlled manner. Spin networks have been proposed as 
quantum channels |21j where the collective dynamics of their low lying excitati- 
ons is exploited for transporting quantum information. But tools from quantum 
information theory also start influencing numerical methods as the density ma- 
trix renormalization group and the design of new efficient simulation strategies 
for many-body systems (see for example |134| 11351 1132j ). Of particular interest 
for this thesis will be the extensive analysis of entanglement in quantum critical 
models [501 1551 ITTT]. 

One important conclusion from the enormous bulk of work concerned with 
entanglement at quantum phase transitions is that multipartite quantum corre- 
lations are typically playing a dominant role [1] . This establishes an interconnec- 
tion with the field of invariant theory, where the quantification and classification 
of multipartite entanglement provides with a plethora of open problems also in- 
teresting in mathematics. Without a clear perspective of which multipartite 
quantum correlations might have relevance for certain complex phenomena in 
condensed matter physics and - most importantly - in absence of a full classifi- 
cation of entanglement, essentially those measures are investigated that can be 
easily computed. Though important in its own right, this pragmatic approach 
waives the main scope behind such an analysis, namely the understanding of the 
underlying entanglement pattern and its interconnection with complex physical 
phenomena. 

Most of this thesis represents a summary of a selection of work published 
during the last few years f90lEll5ll9lll93llSll92lll5llT.lt will start with an 
overview over largely employed entanglement measures followed by a selection of 
results of their analysis for condensed matter systems emphasizing the relevan- 
ce of multipartite entanglement. Thereafter, the concept of genuine multipartite 
entanglement is introduced and an approach for the construction of genuine mul- 
tipartite entanglement measures is presented. Relevant new features appearing 
in the multipartite case, as compared to bipartite measures, are highlighted, 
together with an outline of how to measure multipartite entanglement in pure 
states in the laboratory. 



12 



Kapitel 2 



Pairwise and Bipartite 
Entanglement 



The problem of measuring entanglement is a vast and lively field of research in 
its own. Numerous different methods have been proposed for its quantification 
In this Section we do not attempt to give an exhaustive review of the field. 
Rather do we want to introduce those measures that are largely being used 
to quantify entanglement in many-body systems. Comprehensive overviews of 
entanglement measures can be found in [211 [HI HH Hffl IM [Ml [M|- Also a 
method for detecting entanglement is outlined that is based on entanglement 
witnesses. 

2.1 Bipartite entanglement in pure states 

Bipartite entanglement of pure states is conceptually well understood, although 
quantifying it for local dimensions higher than two still bears theoretical chal- 
lenges |1391 [S3] . A pure bipartite state is not entangled if and only if it can be 
written as a tensor product of pure states of the parts. It is an important fact 
with this respect that for every pure bipartite state | ipAs) (with the two parts, 
A and B), two orthonormal bases {| ipA,i)} and {\4>B,j)} exist such that | ipAs) 
can be written as 



where are positive coefficients. This decomposition is called the Schmidt 
decomposition and the particular basis coincide with the eigenbasis of the cor- 
responding reduced density operators 



The density operators pA and ps have common spectrum, in particular are they 
equally mixed. Since only product states lead to pure reduced density matrices, 
a measure for their mixedness points a way towards quantifying entanglement in 




(2.1) 
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this case. Given the state | iPab), we can thus take its Schmidt decomposition, 
Ea. (|2.ip . and use a suitable function of the ai to quantify the entanglement. 

It is interesting that an entanglement measure E is fixed uniquely after 
imposing the following conditions 

1. i? is invariant under local unitary operations E is indeed a function 
of the aiS only). 

2. i? is continuous (in a certain sense also in the asymptotic limit of infinite 
copies of the state; see e.g. Ref. |101j ). 

3. E is additive, when several copies of the system are present: 

E{\ ^Pab) ® I 4'ab)) = E{\ i,AB)) + E{\ cpab))- 

The unique measure of entanglement satisfying all the above conditions is the 
von Neumann entropy of the reduced density matrices 

S{pa) - S{pb) = - E ' (2.2) 

■i 

this is just the Shannon entropy of the moduli squared of the Schmidt coeffi- 
cients. In other words: one possible answer on the question of how entangled a 
bipartite pure state is, can be given by the von Neumann entropy of (either of) 
the reduced density matrices. The amount of entanglement is generally difficult 
to define once we are away from bipartite states, but in several cases we can 
still gain some insight into many-party entanglement if one considers different 
bipartitions of a multipartite system. In particular if no reduced density matrix 
is pure, then the state is called globally entangled. 

It is worth to notice that a variety of purity measures are admissible when the 
third condition on additivity is omitted. In principle, there are infinitely many 
measures for the mixedness of a density matrix; two of them will typically lead to 
a different ordering when the Hilbert space of the parts has a dimension larger 
than two. This is essentially equivalent to saying that different inequivalent 
classes of entanglement exist in these cases. In contrast, if we trace out one of 
two qubits the corresponding reduced density matrix pA contains only a single 
independent and unitarily invariant parameter: its smallest eigenvalue. This 
implies that each monotonic function [0,1/2] i— > [0,1] of this eigenvalue can 
be used as an entanglement measure. Though, also here an infinity of different 
mixedness measures exists, here all lead to the same ordering of states with 
respect to their entanglement, and in this sense all are equivalent. A relevant 
example is the (one-) tangle [5Dj 

Ti[pA]:=AdeipA. (2.3) 
By expressing pA in terms of spin form factors 




where {S") = trs^PAS"') and 5" = ^ct", cr" {a = x,y,z} being the PauH 
matrices, it follows that 

Tr[pA] = l-^{{S-f + {Syf + {S^f). 
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(2.5) 



applies, where C is called the pure state concurrence [59 lll45j . and * indicates the 
complex conjugation in the eigenbasis ola^ . It is worth emphasizing already here 
that all measures of pairwise qubit entanglement can hence be expressed in terms 
of the modulus squared of the expectation value of an antilinear operator. This 
innocent looking detail enhances the minimally required invariance with respect 
to local SU (2) transformations to the local invariance group SL(2). The latter is 
known to be the relevant invariance group for the classification of entanglement 
into SLOCC classes, where generalized local measurements are admitted in a 
probabilistic way j38| . This observation will be a key element, paving the way 
towards the quantification and classification of multipartite entanglement. 
The von Neumann entropy can be expressed as a function of the (one-) tangle 



where h{x) —: — xlogj x ~ {I — x) log2(l — x) is the binary entropy. Both the 
tangle and the concurrence lead to the important monogamy inequality [30l [89] 
which will be discussed in the next Section. 

2.2 Pairwise qubit entanglement in mixed states 

Subsystems of a many-body (pure) state will generally be in a mixed state, 
and then even different concepts of entanglement do exist. Three important re- 
presentatives are the entanglement cost Ec, the distillable entanglement Eu 
(both defined in Ref. [T7]) and the entanglement of formation Ep [THj. Where- 
as Ej} and Ec are asymptotic limits of multi-copy extraction probabilities of 
Bell states and creation from such states, the entanglement of formation is the 
amount of pure state entanglement needed to create a single copy of the mixed 
state. Very recently, a proof of the full additivity of Ep has been presented [96] . 
which impHes that for bipartite systems both concepts coincide (see e.g. |136| ). 
i.e. En = Ec- The conceptual difficulty behind the calculation of Ep lies in 
the infinite number of possible decompositions of a density matrix. Therefore, 
even knowing how to quantify bipartite entanglement in pure states, we cannot 
simply apply this knowledge to mixed states in terms of an average over the 
mixtures of pure state entanglement. The problem is that two decompositions 
of the same density matrix usually lead to a different average entanglement. 
Which one do we choose? It turns out that we must take the minimum over 
all possible decompositions, simply because if there is a decomposition where 
the average is zero, then this state can be created locally without need of any 
entangled pure state, and therefore Ep — 0. The same conclusion can be drawn 
from the requirement that entanglement must not increase on average by means 
of local operations including classical communication (LOCC). A minimal set 
of requirements every entanglement measure has to fulfill has lead to the notion 
of an entanglement monotone |133j . 

The entanglement of formation of a state p is therefore defined as 





(2.6) 
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where the minimum is taken over all realizations of the state pab = X]j -Pi I V'i) i^j I 
and S{pA,j) is the von Neumann entropy of the reduced density matrix pAj '■= 
tr B IV'i) (V'i I- Ea. (|2.6p is the so-called convex roof of the entanglement of for- 
mation for pure states, and a decomposition leading to this convex roof value is 
called an optimal decomposition. 

For systems of two qubits, an analytic expression for Ep does exist and it is 
given by 

SH.)^-i:^^±|^.nV^±|^ (2.7, 

tT=± 

where C (p) is the convex roof of the pure state concurrence [145|, I146j which 
has been defined in the previous section. Its convex roof extension is encoded 
in the positive Hermitean matrix 

R = ^pp^p = ^p{ay ® ay)p* [ay ® ay)^p , (2.8) 

with eigenvalues > • • • > A4 in the following waj0 

C = max{Ai - A2 - A3 - A4, 0} . (2.9) 

As the entanglement of formation is a monotonous function of the concurrence, 
also C itself or its square T2 can be used as entanglement measures. This is 
possible due to a curious peculiarity of two-qubit systems: namely that a con- 
tinuous variety of optimal decompositions exist [145j . The concurrence C and 
the tangle ti both range from (no entanglement) to 1. 

By virtue of (|2.8p and (|2.9p . the concurrence in a spin-1/2 chain can be 
computed in terms of up to two-point spin correlation functions. For simplicity 
we consider a case where the model has a parity symmetry [1. For this case the 
reduced density matrix pij for spins placed at sites i and j assumes a simple 
form in the basis {| 00) , | 01) , 1 10) , | 11)} 1 



(2) 
Pi3 



I a^j Cij \ 

Xij Zij 

z*, y^j 



(2.10) 

with real a, 6, x, y, and complex c, z. The concurrence results to be 



Cij = 2max{0, |cy | - ^/x^jyij, |%| - \/(Mjb~} . (2-11) 

For translational invariant systems: x = y; for real Hamiltonians and stationary 
states: c, z e IR. Each entry of the matrix pij is then simply related to one- and 
two-point correlation functions, 

= 2 max {O, C/^, C//} . (2.12) 



with its time reversal. 

■^In this case all the components of the wave function have an even (or odd) number of 
flipped spins 

^when discussing qubits or spin-1/2 systems both notations | 0) , | 1) and | T) > I i) will be 
used for the eigenstates of S^. 
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where 



= W + 9!!\-\j(^\+9!^J -M^ , (2.13) 

= W-9m+9tf~\, (2.14) 

with g"" = (5'f S'") and Mz = (S^) (assuming translational invariance). A state 
with dominant fideUty of parallel and anti-parallel Bell states is characterized 
by dominant and C^^ , respectively. This was shown in [46], where the con- 
currence was expressed in terms of the fully entangled fraction as defined in [19j . 

The importance of the tangle and the concurrence is due to the monogamy 
inequality derived in [30] for three qubits. This inequality has been proved to 
hold also for n-qubits system In the case of many-qubits (the tangle may 
depend on the site i) it reads 

J2C-,<ri,^. (2.15) 

The so called residual tangle Ti^i — X^j^z ^fj^ then a measure for multipartite 
entanglement [30, 89J not stored in pairs of qubits only. We finally mention 
that the antilinear form of the concurrence was the key for the first explicit 
construction of a convex roof, and hence its extension to mixed states [SH I145| 

[m]. 

Another measure of entanglement we mention is the relative entropy of 
entanglement |127j . It can be applied to any number of qubits in principle (or 
any dimension of the local Hilbert space). It is formally defined as E{a) :— 
miupgp S{a\ \p), where S(a\ \p) = tr cr [Incj — Inp] is the quantum relative entro- 
py. This relative entropy of entanglement quantifies the entanglement in a by 
its distance from the set V of separable states. The main difficulty in computing 
this measure is to find the disentangled state closest to p. This is in general a 
difficult task, even for two qubits. In the presence of certain symmetries - which 
is the case for e.g. eigenstates of certain models - an analytical access is possible. 
In these cases, the relative entropy of entanglement becomes a very useful tool. 
The relative entropy reduces to the entanglement entropy in the case of pure bi- 
partite states; this also means that its convex roof extension coincides with the 
entanglement of formation, and is readily deduced from the concurrence [145] . 

2.3 Localizable entanglement 

A different approach to entanglement in many-body systems arises from the 
quest to swap or transmute different types of multipartite entanglement into 
pairwise entanglement between two parties by means of generalized measures 
on the rest of the system. In a system of interacting spins on a lattice one could 
then try to maximize the entanglement between two spins (at positions i and j) 
by performing measurements on all the others. The system is then partitioned 
in three regions: the sites i, j and the rest of the lattice. This concentrated 
pairwise entanglement can then be used e.g. for quantum information processing. 
A standard example is that the three qubit GHZ state (1/V2)(| 000) -I- | 111)) 
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after a projective measure in x-direction on one of the sites is transformed into 
a two qubit Bell state. 

The concept of localizable entanglement has been introduced in [1311 1103j . 
It is defined as the maximal amount of entanglement that can be localized, on 
average, by doing local measurements in the rest of the systenfl. In the case of 
N parties, the possible outcomes of the measurements on the remaining N — 2 
particles are pure states lips) with corresponding probabihties . The localizable 
entanglement Eioc on the sites i and j is defined as the maximum of the average 
entanglement over all possible outcome states | '>jjs)ij 

EioS,j) = supf Y.P-^^\ (2.16) 

where £ is the set of all possible outcomes {ps, IV's)) of the measurements, and 
E is the chosen measure of entanglement of a pure state of two qubits (e.g. the 
concurrence). Although very difhcult to compute, lower and upper bounds have 
been found which allow to deduce a number of consequences for this quantity. 

An upper bound to the localizable entanglement is given by the entanglement 
of assistance [HH] obtained from localizable entanglement when also global and 
joint measurements were allowed on the N — 2 spins E|. A lower bound of the 
localizable entanglement comes from the following theorem |131j 

Theorem 2.3.1 Given a (pure or mixed) state of N qubits with reduced cor- 
relations Q°'j^ = (Sf'Sj) — {Sf'){Sj) between the spins i and j and directions 
a and (3 then there always exist directions in which one can measure the other 
spins such that this correlation do not decrease, on average. 

It then follows that a lower bound to localizable entanglement is fixed by the 
maximal correlation function between the two parties (one of the various spin- 
spin correlation functions Q^'^'jfl. 

2.4 Entanglement witnesses 

It is important to realize that not just the quantification of many-party entan- 
glement is a difficult task; it is an open problem to tell in general, whether a 
state of n parties is separable or not although a formal solution of the problem 
can be written in several forms. It is therefore of great value to have a tool 
that is able to merely certify if a certain state is entangled. An entanglement 
witness is a Hermitean operator which is able to detect entanglement in a 
state. The basic idea is that the expectation value of the witness W for the 
state p under consideration exceeds certain bounds only when p is entangled. 
An expectation value of W within this bound however does not guarantee that 
the state is separable. Nonetheless, this is a very appealing method also from an 

^ These operations in principle need not to be local operations in terms of the multipartite 
setting of single spins on all the chain. 

^The entanglement of assistance is defined as the maximum average entanglement among 
the pure state realizations of the state under consideration |17) 

^It has been argued recently |53II52I that in order to extend the entanglement of assistance 
and the localizable entanglement to being an entanglement monotone |133j one should admit 
also local operations including classical communication on the extracted two spins, this was 
named entanglement of collaboration. 
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experimental point of view, since it is sometimes possible to relate the presence 
of the entanglement to the measurement of few observables. 

Simple geometric ideas help to explain the witness operator W at work. Let T 
be the set of all density matrices and let £ and P be the subsets of entangled and 
separable states, respectively. The convexity of 2? is a key property for witnessing 
entanglemenl[3. The entanglement witness is then an operator defining a hyper- 
plane which separates a given entangled state from the set of separable states. 
The main scope of this geometric approach is then to optimize the witness 
operator |71j or to replace the hyper-plane by a curved manifold, tangent to the 
set of separable states [56] H. We have the freedom to choose W such that 

tr (pdW) < 

for all disentangled states pD G 2^- Then, 

tr (pW) > 

implies that p is entangled. A caveat is that the concept of a witness is not 
invariant under local unitary operations (see e.g. [H]). 

Entanglement witnesses are a special case of a more general concept, name- 
ly that of positive maps. These are injective superoperators on the subset of 
positive operators. When we now think of superoperators that act non-trivially 
only on a sub-Hilbert space, then we may ask the question whether a positive 
map on the subspace is also positive when acting on the whole space. Maps that 
remain positive also on the extended space are called completely positive maps. 
The Hermitean time evolution of a density matrix is an example for a comple- 
tely positive map. Positive but not completely positive maps are important for 
entanglement theory [611 IM] 

Theorem 2.4.1 A state pab is entangled if and only if a positive map A exists 
(not completely positive) such that 

{1a<»Ab)pab < . 

For a two dimensional local Hilbert space the situation simplifies considerably 
in that any positive map P can be written as P = CPi + CP2TB , where 
CPi and CP2 are completely positive maps and Tb is a transposition operation 
on subspace B. This decomposition tells that for a system of two qubits the 
lack of complete positivity in a positive map is due to a partial transposition. 
This partial transposition clearly leads to a positive operator if the state is a 
tensor product of the parts. In fact, also the opposite is true: a state of two 
qubits PAB is separable if and only if > that is, its partial transposition 
is positive. This is very simple to test and it is known as the Peres-Horodecki 
criterion [HI] . The properties of entangled states under partial transposition 
lead to a measure of entanglement known as the negativity. The negativity Nab 
of a bipartite state is defined as the absolute value of the sum of the negative 
eigenvalues of Pab- "The logarithmic negativity is then defined as 

En ^\og2 2{2NAB + l)- (2.17) 

^This is based on the Hahn-Banach separation theorem, stating that given a convex set 
and a point outside there exists a plane that separates the point from the set. 
®For other geometric aspects of entanglement see 1651 1161 170) . 
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For bipartite states of two qubits, p/^^ has at most one negative eigenvalue |110j . 
For general multipartite systems and higher local dimension there are entangled 
states with a positive partial transpose, known as bound entangled states [Tll62j. 

2.5 Indistinguishable particles 

There is an ongoing debate on which definition of entanglement for indistin- 
guishable particles will be the most useful from a physical point of view. This 
uncertainty is responsible for the vast variety of quantities studied, when the 
entanglement of itinerant fermion and boson systems is discussed. 

The problem is that for indistinguishable particles the wave function is 
(anti-) symmetrized and therefore the definition of entangled states as given 
in the previous Section does not apply. In particular, it does not make sen- 
se to consider each individual particle as parts of the partition of the system. 
Following |491 147j one can address the problem of defining entanglement in an 
ensemble of indistinguishable particles by seeing if one can attribute to each of 
the subsystems a complete set of measurable properties, e.g. momenta for free 
pointless particles. Quantum states satisfying the above requirement are pre- 
cisely the (anti-) symmetrization of a product state of (fermions) bosons, and 
represent the separable states for indistinguishable particles. 

There is another crucial difference between the entanglement of (indistin- 
guishable) spin-1/2 particles and that of qubits. Let us consider two fermions 
on two sites. Whereas the Hilbert space Tig of a two-site spin lattice has di- 
mension AivnTLs = 4, the Hilbert space Tif for two fermions on the same lattice 
has dimension dvcaTLf = 6. This is due to the possibility that both fermions, 
with opposite spins, can be located at the same lattice site. When choosing the 
following numbering of the states 
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(2.18) 



and the definition — fjfj \ 0), there are Bell states analogous to those oc- 
curring for distinguishable particles (| 1, 3) ± | 2, 4))/V2 and (| 1, 4) ± | 2, 3))/V2. 
There are however new entangled states, as (|1,2) ± |3,4))/V2, where both 
fermions take the same position. The local Hilbert space is made of four states 
labelled by the occupation number and the spin, if singly occupied. The site- 
entanglement of indistinguishable particles is then defined as the entanglement 
of the corresponding Fock states. It can be measured e.g. by the local von Neu- 
mann entropy. This quantity is the analogue to the one-tangle for qubits, but 
the local Hilbert space dimension is 4 due to the possibility of having empty 
and doubly occupied sites. Also the quantum mutual information [Hj can be 
defined in this way, quantifying the total amount (classical and quantum) of 
correlations stored in a given state of a second quantized system. 

For spinless fermions a one-to-one mapping to spin-1/2 chains exists in 
one spatial dimension - the Jordan- Wigner transformation. Due to the non- 
locality of the Jordan Wigner transformation, quantitative deviations between 
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the entanglement of spinless fermions and the Jordan- Wigner relative may oc- 
cur for distant sites; as far as nearest-neighbor entanglement is considered, both 
concepts completely coincide. 

Although it is known how the entanglement of indistinguishable particles 
can be quantified, as will be seen in the following part, the major part of the li- 
terature on second quantized systems considers the site-entanglement described 
above or the entanglement of degrees of freedom, singled out from a suitable 
set of local quantum numbers (e.g. the spin of the particle at site i). In both 
cases, entanglement measures for distinguishable particles (see Sections 13.4.11 
and 13.4.3]) can be used. 



2.5.1 Two Fermion entanglement 

In this paragraph we summarize the distinct features appearing in the quantifica- 
tion and classification of Fermion entanglement. Due to the antisymmetry under 
particle exchange, there is no Schmidt decomposition for Fermions. Neverthe- 
less, a Fermionic analogue to the Schmidt rank, which classifies entanglement 
in bipartite systems of distinguishable particles does exist: the so called Slater 
rank. A generic state of two-electrons on two lattice sites can be written as 

4 

where is a 4 x 4 matrix which can be assumed antisymmetric and normalized 
as tr uj^uj = i (or equivalently tr uj*uj ~ ^^)- Since here the local entities whose 
entanglement shall be studied, are the particles, unitary transformations act on 
the 4-dimensional single particle Hilbert space. Due to the indistinguishability of 
the particles, the transformation must be the same for each of the particles. Gi- 
ven a unitary transformation U E SU(4) such that fj := Ujkff^, the transformed 
state is given by | w') where cu' := UojU'^ . The above unitary transformation 
preserves the antisymmetry, and every pure state oj of two spin-1/2 particles on 
two sites can be transformed into the normal form 



(2.20) 



In fact, every two-particle state within a _D-dimcnsional single particle Hilbert 
space can be transformed into the normal form 

= diag{Zi,...,Z„Zo} (2.21) 

- (A o) (2.22) 
(Zo)., = ; 1,3 E{l,...,D~2r} (2.23) 

where r is then called the Slater rank of the pure Fermion state |111|, 11121 HO] ■ 
Following the definition in the introduction to this Section, a pure Fermion state 
is entangled if and only if its Slater rank is larger than 1. 
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It is important to notice that the above concept of entanglement only de- 
pends on the dimension of the Hilbert space accessible to each of the particles 
(this includes indistinguishable particles on a single D- level system). 

For electrons on an L-site lattice the "local" Hilbert space dimension is 2L, 
and the question, whether a pure state living in a 2L-dimensional single particle 
Hilbert space has full Slater rank, can be answered by considering the Pfaffian 

of w [251 [55] 



Pf H := sign(7r) Y[ ^7r(2j-i) 



(2.24) 



which is non-zero only if uj has full Slater rank L. In the above definition 5^ 
denotes those elements tt of the symmetric group S2L with ordered pairs, i.e. 
7r(2m — 1) < 7r(2m) for all m < L and 7r(2fc — 1) < 7r(2m — 1) for k < m. Notice 
that relaxing the restriction to 5^ just leads to a combinatorial factor of 2^L! 
by virtue of the antisymmetry of uj and hence can we write 



pf[u;] 



2L 



(2.25) 



where e-' 



is the fully antisymmetric tensor with 



,2L 



1. There is a 

simple relation between the Pfafhan and the determinant of an antisymmetric 
even-dimensional matrix: pf[ci^]^ = det[Lj]. 

For the simplest case of two spin-1/2 Fermions on two lattice sites the Pfaffi- 
an reads pf [w] = ti^i,2'^3,4— t^i,3'^2,4+t^i,4'^2,3- Normalized in order to range in the 
interval [0, 1] this has been called the Fermionic concurrence C[| w)] [111|,I112| |40] 
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(2.26) 



where 



UJ 



-e'^'^^uj 



k,l 



(2.27) 



has been called the dual to uj. Then, \ uj) =: I? | w) is the analogue to the conjuga- 
ted state in [59l 11451 If 23j leading to the concurrence for qubits. It is important 
to notice that the Pfaffian in Eq. (j2.24| l is invariant under the complcxification 
of su(2L), since it is the expectation value of an antilinear operator, namely the 
conjugation V for the state |a;). Since this invariant is a bilinear expression in 
the state coefficients, its convex roof is readily obtained jl23j by means of the 
positive eigenvalues of the 6x6 matrix 



R = ^pVpV^p . 

The conjugation 2?, expressed in the basis {| 1,2) 
(see Ea. (|2.18p ). takes the form 



(2.28) 
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where C is the complex conjugation. Notice that the center part of this matrix 
(in bold face) is precisely ay (g) ay and indeed corresponds to the Hilbert space of 
two qubits. The remaining part of the Hilbert space gives rise to an entanglement 
of different values for the occupation number. This type of entanglement was 
considered practically useless and has been referred to as the fluffy bunny jl421 
1128] in the literature. 

For a single particle Hilbert space with dimension larger than 4 one en- 
counters similar complications as for two distinguishable particles on a bipartite 
lattice and local Hilbert space dimension larger than 2, i.e. for two qudits. This 
is because different classes of entanglement occur, which are characterized by 
different Slater rank as opposed to their classification by different Schmidt rank 
for distinguishable particles. The Slater rank can be obtained by looking at Pfaf- 
fian minors [85] : if the Slater rank is r, all Pfaffian minors of dimension larger 
than 2r are identically zero. 



2.5.2 Multipartite Entanglement for Fermions 

For indistinguishable particles the only classification available up to now is to 
check whether or not a pure state has Slater rank one. Eckert et al. formulated 
two recursive lemmata |40j 

Lemma 2.5.1 A pure M-Fermion state has Slater rank one if and only if 
has Slater rank one or zero for all a G C^" . 

Lemma 2.5.2 A pure M-Fermion state has Slater rank one if and only if 



2n 



E 



1 M-2 

^il.....,iM-l^jl,---,jM-l'^jl ■ ■ ■ 



-lM-ltM3M-ljMai...a2n-2] 



(2.30) 



for all a\ . . . ,a^-^ S C^" and all < ai < ■ ■ ■ < a2n-2 < 2n. 

They can be summarized as follows: let an A^-clcctron state be contracted with 
N — 2 arbitrary single electron states encoded in the vectors a^ as a|./^ | 0) 
(j — 1, . . . , N ~2 and sum convention) to a two-electron state. Then the Pfafhan 
of the two-electron state is zero if and only if the original state (and hence all 
intermediate states in a successive contraction) has Slater rank one. This means 
that all 4-dimensional Pfafhan minors of lo are zero. 

Instead of the Pfaffian of uj, also the single-particle reduced density matrix 
can be considered, and its von Neumann entropy as a measure for the quantum 
entanglement has been analyzed in [72l [95] . It is important to remind that for 
distinguishable particles the local reduced density matrix has rank one if and 
only if the original state were a product. This is no longer true for indistinguis- 
hable particles. For an A^-particle pure state with Slater rank one the rank of the 
single-particle reduced density matrix coincides with the number of particles, 
N. A measure of entanglement is then obtained only after subtraction of the 
constant value of the von Neumann entropy of a disentangled state. This must 
be taken into account also for the extension of the measure to mixed states. 
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2.5.3 "Entanglement of particles 



Entanglement in the presence of superselection rules (SSR) induced by particle 
conservation has been discussed in Refs. [HI 11431 11141 1113j . The main diffe- 
rence in the concept of entanglement of particles [143j from the entanglement 
of indistinguishable particles as described in the preceding section (but also to 
that obtained from the reduced density matrix of e.g. spin degrees of freedom of 
indistinguishable particles) consists in the projection of the Hilbert space onto a 
subspace of fixed particle numbers for either part of a bipartition of the system. 
The bipartition is typically chosen to be space-like, as motivated from experi- 
mentalists or detectors sitting at distinct positions. E.g. two experimentalists, in 
order to detect the entanglement between two indistinguishable particles, must 
have one particle each in their laboratory. 

This difference induced by particle number superselection is very subtle and 
shows up if multiple occupancies occur at single sites for Fermions with some 
inner degrees of freedom, as the spin. Their contribution is finite for finite dis- 
crete lattices and will generally scale to zero in the thermodynamic limit with 
vanishing lattice spacing. Therefore both concepts of spin entanglement of two 
distant particles coincide in this limit. Significant differences are to be expected 
only for finite non-dilute systems. It must be noted that the same restricti- 
ons imposed by SSR which change considerably the concept of entanglement 
quantitatively and qualitatively, on the other hand enable otherwise impossible 
protocols of quantum information processing jll41 1113j . 

Wiseman and Vaccaro project an iV-particle state \iPn) onto all possible 
subspaces, where the two parties have a well defined number {ua, riB = N — ua) 
of particles in their laboratory [143j . Let be the respective projection, 

and let be the weight ( ^/'[n^] | tp[nA] )/( V'w I ) of this projection. Then 
the entanglement of particles Ep is defined as 

Ep[\ ^n)] = ^P„£;Af[V'M] (2.31) 

n 

where Em is some measure of entanglement for distinguishable particles. Alt- 
hough this certainly represents a definition of entanglement appealing for ex- 
perimental issues, it is sensitive only to situations, where e.g. the two initially 
indistinguishable particles eventually are separated and can be examined one- 
by-one by Alice and Bob. Consequently, "local operations" have been defined 
in |143| as those performed by Alice and Bob in their laboratory after having 
measured the number of particles. 

Verstraete and Cirac pointed out that the presence of SSR gives rise to a 
new resource which has to be quantified. They have proposed to replace the 
quantity Ep with the SSR- entanglement of formation. This is defined as 

Ei^^^\\ipN)] = min V'p„£^M[V'n] 

n 

where the minimization is performed over all those decomposition of the density 
matrix where the | ■!/))„ are eigenstates of the total number of particles |114[I113] . 

2.5.4 Entanglement for Bosons 

The quantification and classification of boson entanglement is very close in spi- 
rit to that of Fermions as described in Section 12.5.11 We therefore will only 
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emphasize the marking differences for bosonic entanglement. 

In the bosonic case the matrix w in Eq. (|2.19p is symmetric under permutati- 
ons of the particle numbers. Consequently, for any two-particle state of indistin- 
guishable bosons, Lo can be diagonalized by means of unitary transformations 
of the single particle basis. This leads to the Schmidt decomposition for bo- 
sons |40) . An curious feature distinguishing this case from the entanglement 
measures of distinguishable particles is that the Schmidt decomposition is not 
unique. In fact, any two equal Schmidt coefficients admit for a unitary trans- 
formation of the two corresponding basis states, such that the superposition of 
the two doubly occupied states can be written as a symmetrized state of two 
orthogonal states [72j [483 • This is the reason why it is not directly the Schmidt 
rank, but rather the reduced Schmidt rank - obtained after having removed all 
double degeneracies of the Schmidt decomposition - that determines whether 
or not a state is entangled. This non-uniqueness of the Schmidt rank is also 
responsible for the ambiguity of the von Neumann entropy or other purity mea- 
sures of the single particle reduced density matrix as an entanglement measure 
for Bosons [15]. 

With Zi being the Schmidt coefficients with degeneracy gi, the reduced 
Schmidt rank is at most ^ + 2 { ^ } , where {.} denotes the non- integer part. As 
a consequence, a Schmidt rank larger than two implies the presence of entangle- 
ment. Schmidt rank 2 with degenerate Schmidt coefficients can be written as a 
symmetrized product of orthogonal states and consequently is disentangled [35] ■ 
This feature is also present in the iV-boson case, where in presence of up to A^- 
fold degenerate Schmidt coefficients the corresponding state can be rewritten as 
a symmetrization of a product. 

For bipartite systems lo has full Schmidt rank if det w ^ 0. A Schmidt rank 1 
can be verified by the same contraction technique described for the Fermion case 
in the previous section, where the Pfaffian must be replaced by the determinant. 
This applies to both the bipartite and the multipartite case [ID]. 
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Kapitel 3 

Entanglement in condensed 
matter systems 



Traditionally many-body systems have been studied by looking for example at 
their response to external perturbations, various order parameters and exci- 
tation spectrum. The study of the ground state of many-body systems with 
methods developed in quantum information might unveil new insight. In this 
Section we classify the properties of the ground state of a many-body system ac- 
cording to its entanglement. We first concentrate on spin systems. Spin variables 
constitute a good example of distinguishable objects, for which the problem of 
entanglement quantification is most developed. Various aspects mainly of pair- 
wise entanglement will be discussed with some short reference on the properties 
of bipartite entanglement - as the block entropy - and a comment on the locali- 
zable entanglement. Then will we change focus onto itinerant fermion systems. 



3.1 Model systems 

The model Hamiltonian for a set of localized spins interacting via nearest neigh- 
bor exchange coupling in a d-dimensional lattice can be written as 



h(7,a,/j7j) = jJ2 



2 



(3.1) 



In the previous expression i,j are lattice points, (•) constraints the sum over 
nearest neighbors and S" {a = x,y,z) are spin-1/2 operators. The nomencla- 
ture of the various model deriving from Eq. (|3.ip is shown in table 13. ll A positive 
(negative) exchange coupling J favors anti-ferromagnetic (ferromagnetic) orde- 
ring in the xy— plane. The parameters 7 and A account for the anisotropy in 
the exchange coupling, h is the transverse magnetic field. 

The ground state of Eq. ()3.1|) is in general entangled, but for any value of the 
coupling constants 7 and A, J > a value h{ for the magnetic field exists in 
d— 1,2, where the ground state is factorized [551ll07j . The so called factorizing 
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XX 








XY 


^0 
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XXZ 





^0 


XYZ 
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Ising 


1 






Tabelle 3.1: Nomenclature of the various models deriving from Eq. (|3.ip 
field hi is given by 

h, = ^Jv/(l + A)2-(^/2)2 

where z is the coordination number of the lattice. Note that the result for the 
factorizing field is rigorous irrespective the integrability of the Hamiltonian. 

In d = 1 the model is exactly solvable in several important cases. In the 
next two paragraphs we illustrate some of the results obtained for the exactly 
solvable transverse XY model (A = and < 7 < 1). The quantum Ising 
model is a special case corresponding to 7 = 1 while the (isotropic) XX-moAel 
is obtained for 7 = 0. In the isotropic case the model possesses an additional 
symmetry, resulting in the conservation of the magnetization along the z-axis. 
For any value of the anisotropy the model can be solved exactly [7S1 [Ml US] 
by a Jordan- Wigner- and a successive Bogoliubov transformation. The Lipkin- 
Meshkov-Glick model |74| which also will appear in the discussion, emerges from 
the transverse XY-models when the range of spin-exchange, or equivalently the 
connectivity, is set to infinity. 

The properties of the Hamiltonian are governed by the dimensionless coup- 
ling constant A = J/2h. In the interval < 7 < 1, the system undergoes a 
second order quantum phase transition at the critical value Ac = 1 [109| I122j . 
The order parameter is the in-plane magnetization (e.g. in a;-direction: (5*^)), 
which is different from zero for A > 1. The magnetization along the z-direction, 
{S^), is different from zero for any value of A with singular behavior of its first 
derivative at the transition. This is reflected also in the singularity present in 
the second derivative of the ground state energy with respect to A. In the whole 
interval < 7 < 1 the transition belongs to the Ising universality class. For 
7 = the quantum phase transition is of the Berezinskii-Kosterlitz-Thouless 
type. 

3.2 Bipartite entanglement and quantum phase 
transitions 

Many scientific investigations have been devoted to the study of entanglement 
close to quantum phase transition (QPT). In contrast to a standard thermody- 
namic phase transition, a QPT is a phenomenon that occurs at zero temperature. 
Its essence consists in a significant qualitative change of the ground state of a 
model Hamiltonian, which is induced by the change of an external parameter or 
coupling constant [109j . The main idea is that this drastic change of the ground 
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state should be accompanied by characteristic entanglement patterns. Similar 
to standard phase transition, also a quantum critical point is characterized by 
a diverging correlation length ^, which is responsible for the singular behavior 
of different physical observables. The critical properties in the entanglement we 
are going to summarize below admit a screening of the qualitative change of 
the state of the system experiencing a quantum phase transition. In order to 
avoid possible confusion, it is worth to stress that the study of entanglement 
close to quantum critical points is not supposed to provide new insight into 
the scaling theory of quantum phase transitions; rather, it may be useful for 
a deeper characterization of the ground state wave function of the many-body 
system undergoing a phase change. In the following, we shine some light on the 
behavior of the pairwise entanglement, with a subsequent glance at the block 
entropy 0- 

Pairwise entanglement close to quantum phase transitions was originally 
analyzed in [SHI EH] for the quantum XY model in transverse magnetic field 
in one spatial dimension. For the quantum Ising chain, the concurrence tends 
to zero for A ^ 1 and A <C 1, where the ground state of the system is fully 
polarized along the a;-axis and the z-axis, respectively. Whereas the full pola- 
rization in z-direction for small A is guaranteed by the large magnetic field in 
z-direction, this is not the case in the opposite limit for generic values of 7, and 
the polarization in this case is an effect due to symmetry breaking. A particu- 
larly surprising observation is the short range of the concurrence, in particular 
at the critical point, notwithstanding the diverging range of two-point spin cor- 
relations: the concurrence is zero unless the two sites are at most next-nearest 
neighbors. This short range in the pairwise entanglement is observed in many 
different models [120U1211 l4]. also in higher spatial dimensions. It hence seems to 
be a generic feature rather than a curious exception. This indicates that pair- 
wise entanglement typically plays a secondary role in Cavour of multipartite 
entanglement, as far as quantum phase transitions are concerned. 

In the Ising case, the concurrence is a smooth function of the coupling with its 
maximum well separated from the critical point (see the right inset of Fig l3.ip . 
In contrast, the convex roof of the von Neumann entropy shows a pronounced 
cusp at the critical point j88j. The critical properties of the ground state are 
instead well reflected in the derivatives of the concurrence as a function of A. 
The results for systems of different size (including the thermodynamic limit) 
are shown in Fig l3.1l For the infinite chain ^^^(l) diverges on approaching the 
critical value as 

9AC(l)~^ln|A-A,| . (3.2) 

For finite system size, the precursors of the critical behavior can be analyzed 
by means of finite size scaling. In agreement with the scaling hypothesis, the 
concurrence depends only on the combination N^/'^ {\—\m) in the critical region, 
with critical exponent v = \\ X„i is here the position of the minimum (see the 
left inset of Fig |3.ip . In the case of log divergence the scaling ansatz has to 
be adapted for taking care of the critical characteristic log divergence in the 

^ QPTs were also studied by looking at quantum fidelity |32l I149| or the effect of single bit 
operations |51l 150) 
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Abbildung 3.1: The change in the ground state wave- function in the critical 
region is analyzed considering the derivative of the nearest neighbor concurrence 
as a function of the reduced coupling strength. Different curves correspond to 
different lattice sizes. On increasing the system size, the minimum gets more 
pronounced. Also the position of the minimum changes and tends as N~^'^^ 
(see the left side inset) towards the critical point where for an infinite system a 
logarithmic divergence is present. The right hand side inset shows the behavior 
of the concurrence itself for an infinite system. The maximum is not related to 
the critical properties of the Ising model. [From [50] ] 



quantum Ising universality class 

dxC{l){N,X) - 9aC(1)(7V,Ao) 

- Q[N^/''d,n{X)]~Q[N^^''S,n{Xo)] (3.3) 

where Aq is some non critical value, (5m (A) = A — A,„ and Q{x) ^ Q{oo) lux (for 
large x). Similar results have been obtained for the XY models in this univer- 
sality class (SOj . Although the concurrence describes short-range properties, the 
typical scaling behavior for continuous phase transitions emerges. The analysis 
of the finite size scaling in the so called period-2 and period-3 chains, where the 
exchange coupling varies every second and third lattice sites respectively, leads 
to the same scaling laws in the concurrence |150) . 

Spontaneous symmetry breaking can influence the entanglement in the ground 
state. To see this, it is convenient to introduce the thermal ground state in the 
limit T ^ 

Po^l {\9s"){gs"\ + \9s^){9s^\) = I {\9s-){9s-\ + \gs+){gs+\) . (3.4) 

The symmetry broken states gs~^ and gs~ , which give the correct order parame- 
ter of the model, are superpositions of the degenerate parity eigenstates gs° and 
gs'^. This is the essence of the symmetry breaking in the transverse XY models. 
Being convex, the concurrence in gs^ will be larger than for gs°/'^ [92 . The 
opposite is true for the entropy of entanglement (see Ref. [88^ for the single spin 
von Neumann entropy). It was shown, however, that the effect of the sponta- 
neous parity symmetry breaking does not affect the concurrence in the ground 
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state if it coincides with , Ea. (|2.13p : that is, if the spins are entangled in 
an anti- ferromagnetic way [120j . For the quantum Ising model, the concurrence 
coincides with for all values of the magnetic field, and therefore, the con- 
currence is unaffected by the symmetry breaking, the hallmark of the present 
QPT. For generic anisotropies 7 instead, also the parallel entanglement C^^ is 
observed precisely for magnetic fields larger than the factorizing field [HI]; this 
interval includes the critical point 0. An interpretation of this is that for the 
Ising universality class, the concurrence is insensitive to the symmetry breaking 
close to the critical point, and hence won't play a relevant role in driving this 
transition. This changes at 7 = 0, where the concurrence indeed shows an in- 
finite range. Below the critical field, the concurrence is enhanced by the parity 
symmetry breaking [92^ , as shown in Fig. 13.21 




h/h, ' ' h/h, 



Abbildung 3.2: Left panel: The nearest neighbor concurrence for a chain of 199 
sites is shown for three different values 7 — 0.3, 0.5 and 0.7 (circles, squares 
and diamonds). Full symbols give the results for the even parity ground state. 
Right panel: Difference of the n.n. concurrence with and without broken parity 
symmetry as a function of the transverse field. The maximum relative deviation 
amounts to around 10%; for sufficiently small h, it decreases with 7. Inset: Finite 
size scaling for 7 = 0.7 and h/hf = 0.8 (diamonds) and limiting value (dashed 
line). [From ^] 

Several works were devoted to an entanglement analysis close to this fac- 
torizing field. The point at which the state of the system becomes separable 
comes with an exchange of parallel and anti-parallel sector in the ground state 
concurrence (see Eas. (|2.13p and (|2.14p ). Furthermore, it is observed that the 
range of the concurrence diverges close to the factorizing field, in contrast to 
the typically encountered short range of the concurrence at the critical point. 
There, the range R is taken as the distance of two qubits, beyond which the 
concurrence is zero. For the XY model the range was found to diverge close to 
the factorizing field as [3] 

R(x(^n^^^ \n\X-^ - XJ^\-^ (3.5) 

The divergence of R suggests, as a consequence of the monogamy of the entan- 
glement [301 [M] , that the role of pairwise entanglement is enhanced while ap- 
proaching the separable point [105^ I106[ 1107) . Indeed, for the Ising model (i.e. 

^To see this, it is enough to reaUze that a product state of single spin states is a parity 
eigenstate only if all spins point up or down. 
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7 = 1), the ratio T2/ri — *■ 1 [3], when the magnetic field approaches the factori- 
zing field h{. For 7 7^ 1 and hf < < he it was found that T2/T1 monotonically 
increases with hf and that (T2/ri)|^+ increases with 7^1. The exi- 

stence of a factorizing field emerged as a generic feature of spin chains both for 
short 11051 1106j and long ranged interactions [SH] ; in all these cases the range 
of the two-site entanglement was observed to diverge. 

The section shall be closed with a brief remark concerning the localizable 
entanglement. From its definition and the fact that it is bound from below by 
the largest two-point correlation function, it has been presented as a curiosity 
that an infinite range of entanglement can be present although the correlation 
functions have finite range. As an example for such behavior the AKLT model, 
a spin-1 model which exhibits a Haldane gap, has been presented [103j . The 
infinite range of the localizable entanglement means that by virtue of suitable 
local transformations, different classes of entanglement can be accumulated on 
two sites (see also |102j in this context). At the same time, however, also the 
classical correlations are accumulated and hence would also these localizable 
classical correlations have an infinite range. The curiosity of having infinite range 
entanglement but only finite range correlations is hence only due to looking at 
qualitatively different and not comparable quantities. 

3.3 Dynamics of entanglement 

The interest in studying the properties of entanglement in many-body systems 
has been directed also to the understanding of its dynamical behavior. Entangle- 
ment dynamics has been studied from different perspectives. In a spirit similar 
to the study of propagation of excitations in condensed matter systems, se- 
veral works analyzed the propagation of entanglement starting from a given 
initial state where the entanglement has been created in a given portion of the 
many-body system. One can imagine for example to initialize a spin chain such 
that all the spins are pointing upwards except for two neighbor spins which are 
entangled. Due to the exchange interaction, the initially localized entanglement 
will spread. This propagation may be ballistic in clean systems or diffusive if 
some weak disorder is present. Entanglement localization and chaotic behavior 
could also be observed. An alternative approach is to start with the ground state 
of a Hamiltonian Hq and then let the Hamiltonian change in time. 

Since we are dealing with interacting systems, entanglement can be genera- 
ted or it can change its characteristics during the dynamical evolution. Besides 
the interest in their own, attention to these questions has been also motivated 
by the potential use of one-dimensional spin systems as quantum channels [21j . 
Another important aspect of entanglement dynamics is the possibility to gene- 
rate entangled states with certain desired properties by those interactions that 
are naturally present in a many-body system. This generalizes the setup where 
a Bell state is created by letting two qubits interact for a fixed time by means of 
an exchange coupling of the XX type. In the same spirit one can think of gene- 
rating three-bit entangled GHZ or W states | VF) ~ 1 100) + \ 010) + \ 001) - but 
also other multipartite entangled states - by tailoring the exchange couplings 
in spin networks or their quantum optical counterpart of atoms in an optical 
trap potential. Cluster states are a prominent example for genuinely multipar- 
tite entangled states, which are generated by Ising type two-spin interactions 
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from a fully polarized initial state |58j . 

The most simple situation, which we consider first, is the propagation of 
entanglement in the one-dimensional XX-model, i.e. 7 = and A = in 
Eq. p.ip [St lll9j . This corresponds to free spinless fermions on a discrete lattice. 
Suppose that the initial state of the chain is 

|vl/±(i = 0))^i=«±aj)|0,...0), (3.6) 

namely all the spin are in a fully polarized state except the two at positions i 
and j, which are prepared in one of the Bell states \^±) = 2~^/^{\01) ± |10)). In 
this case the problem has a simple analytical solution. The total magnetization 
is conserved, and the evolution is confined to the sector with only a single spin 
pointing up (single-magnon states). The state of a periodic chain at later times 
is 




(i—l) ^iiJt cos 



ith 



In the thermodynamic limit, TV — > oo, the coefficients can be expressed in terms 
of Bessel functions J„ (x) 

w^fit) = {j.-iim ± (-^)^■-V,_K4J^)} . (3.9) 

The concurrence between two sites located at positions n and m (initially, only 
the sites i and j were maximally entangled) is given as 



Q^(±,t) = 2u;i^f (t)u;2;^';(i) . (3.10) 



and is shown in Fig l3.3l for sites which are symmetrically arranged around 
the initial position of the Bell state \tp±). The Hamiltonian time evolution of 
the XX model leads to a propagation of the single flipped spin through the 
chain. The speed of the propagation is the spin wave velocity in the chain. The 
information exchange or entanglement propagation over a distance of d lattice 
constants approximately takes the time t ~ Hd/ J. The result is independent of 
the external field h^, since the magnetization in 2;-direction is a constant of the 
motion. 



This entanglement wave is also the main feature in the behavior of the von 
Neumann entropy Sn}m of the two sites (n, m) (see Fig l3.4p . 

Interesting additional features appear in the quantum XY model, i.e. for 
7 ^ 0. In this case the magnetization is no longer a constant of the motion (two 
spins can be flipped simultaneously). The calculations were done analytically [S] 
resorting on the exact calculation of correlation functions out of equilibrium [S] . 
The most notable difference in the two-site entanglement is an entanglement 
production from the fully polarized vacuum state. This occurs uniformly along 
the chain and is superposed onto the entanglement wave discussed before (see 
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Abbildung 3.3: Concurrence between sites n = —x,m = x, symmetrically placed 
with respect to the sites i — —I and j ~ 1, where the singlet was initially created. 
[From [ej]. 




Abbildung 3.4: The two-site entropy {t) for pairs of sites symmetrically 

displaced with respect to the initial singlet position (i, j) — (0, 1). The time is 
given in units of the exchange coupling 4 J and the space coordinate x in units 
of the lattice constant. [From [6j]. 
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Abbildung 3.5: The entanglement wave for the nearest neighbor concurrence is 
shown here for the Ising model and A — 0.5. The initial state is a maximal- 
ly entangled state of nearest neighbors on top of the fully polarized vacuum 
state (left) and the ground state (right) respectively. Whereas the propagation 
velocity is unaffected by the initial state, the entanglement wave reduces the 
background entanglement in the ground state of the chain (right). [From [5]] 

left panel of Fig l3.5p . The propagation velocity of the entanglement is almost 
independent on the anisotropy parameter 7 and is in well agreement with the 
sound velocity of the system [118'. What is strongly dependent on 7 is the 
damping coefficient of the entanglement wave: as the anisotropy approaches the 
Ising point 7 = 1, the wave is strongly damped and vanishes already after few 
time units (^ J^^) at the critical coupling. 

The Hamiltonian dynamics primarily generates multipartite entanglement in the 
chain; this can be seen from the residual tangle [3Di ,89^ , which is a measure for 
entanglement of not only pairs of spins. It is seen from Fig. 13.61 that the major 
part of the entanglement in the chain is indeed of multipartite origin. This 
should be expected rather than surprising on the background that multipartite 
entangled cluster states are created from a fully polarized state by means of a 
two-spin Ising- type interactions. 



3.4 Entanglement in second quantized systems 

The theory of entanglement for indistinguishable particles conceptually differs 
from that for distinguishable constituents; the only factoring state of identical 
particles is that of bosons all being in the very same state. 

The major part of the physical applications still concentrates on those mea- 
sures applicable for the entanglement encoded in a certain choice of - distin- 
guishable - quantum numbers (see 12. 5p . This analysis captures, of course, a 
partial aspect of the entanglement encoded in these system; alas, some marking 
peculiarities due to the indistinguishability of the particles are ignored in this 
approach. We will report on the established results from these studies and refer 
to section [2.51 for their relation to fermionic entanglement. 
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Residual tangle: y= A, =1 




Abbildung 3.6: The residual tangle in case of the Ising model and critical coup- 
ling created from the fully polarized vacuum. Multipartite entanglement is clear- 
ly dominant over pairwise entanglement. [From [91j ] 
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3.4.1 Free fermions 



The site-based entanglement of spin degrees of freedom through the Jordan- 
Wigner transformation of spinless fermions has been exploited for calculating the 
concurrence of nearest neighbor sites and the single site von Neumann entropy 
(see Section 12. 5p for the one-dimensional tight-binding model in presence of a 
chemical potential [148] . This model is equivalent to the isotropic XX model in 
a transverse magnetic field. In this specific case, no double occupancy can occur 
and the concept of entanglement coincides with that for spins 1/2. It was found 
that the nearest neighbor concurrence of the ground state at T = is maximal 
at half filling. 

The continuous limit of the tight-binding fermion model is the ideal Fermi 
gas. In this system, the spin entanglement between two distant particles has 
been studied in [124j . There, depending on the dimensionality, the pairwise spin- 
entanglement of two fermions has been found to decrease with their distance 
defining a finite range R of the concurrence. The two spin reduced density matrix 
is 



1 

Pl2 = 



2/2 



1-/^ \ 

1 -p 

0-/2 1 

V 1 -/2/ 



(3.11) 



where f{x) = d^^^-^ with d e {2, 3} being the space dimension and Ji the 
(spherical for d=3) Bessel function of the first kind [1241 l87j. This density matrix 
is entangled if /^ > 1/2. As a consequence, there is spin entanglement for two 
fermions closer than do ~ ^-^^JJ ^^^r d = 3 and do ~ O-^^ij ^'^^ d — 2 {kf is the 
Fermi momentum). A finite temperature tends to diminish slightly the range of 
pairwise spin entanglement [ST]. 

It should not be surprising that non-interacting particles are spin-entangled 
up to some finite distance. It is true that the ground state and even an ar- 
bitrary thermal state of non-interacting fermions has vanishing entanglement 
among the particle^, since the corresponding states are (convex combinations 
of) anti-symmetrized product states. However, disentanglement in momentum 
space typically leads to entanglement in coordinate space. A monochromatic 
plane wave of a single particle for example corresponds to a W state, which 
contains exclusively pairwise entanglement in coordinate space for an arbitrary 
distance of the sites. Furthermore does a momentum cut-off ed kf correspond 
to a length scale of the order kj^. 

It is interesting that a fuzzy detection of the particles in coordinate space 
increases the entanglement detected by the measurement apparatus. In ref. [2^ . 
the two-position reduced density matrix defined by 

P^',u' = {^Ar'y^tiry^Ar')^s{r)) (3.12) 
has been calculated for blurred field operators 

^,{r):^ Jdr'dpil;s{p)D{r-r')e'P''' (3.13) 

^Due to particle- hole symmetry, the concurrence is symmetric around half- filling 
*This statement should not be confused with the non-vanishing entanglement of partic- 
les |143l as observed in |37l . 
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where D{r — r') ~ -^i- exp — ■i^^jS"'' is a- Gaussian distribution describing the 
inaccuracy of the position measurement. This could be understood from the 
blurred field operators being coherent sums of local field operators; the entan- 
glement measured by the apparatus as described above, is the bipartite entan- 
glement between the two regions of width a around r and r' . This entanglement 
is larger than the average of all pairwise contributions out of it due to the super- 
additivity of the entropy/negativity. An analysis in [126j for the three fermion 
spin density matrix revealed that the state carries entanglement within the W- 
class [35] , provided the three particles are in a region with radius of the order of 
the inverse Fermi momentum; a similar reasoning applies to n fermions in such 
a region [T^ITTj. 



3.4.2 su{2) degrees of freedom of interacting fermions 

Itinerant systems, where the focus of interest is the entanglement of degrees 
of freedom forming a representation of su{2) in terms of the fermionic opera- 
tors have been also the subject of investigation. This line has been followed in 
Refs. |148| I117[ I125|. HH I126j for analyzing a connection to BCS superconducti- 
vity and the phenomenon of ry-pairing, a possible scenario for high Tc supercon- 
ductivity. Such states appear as eigenstates of the Hubbard model |147j which 
carry off diagonal long range order. A simplified model of BCS-like pairing for 
spinless fermions has been studied in [148j . The concurrence of the two qubits 
represented by the modes k and —k has been found to be a monotonically in- 
creasing function of the order parameter; it drops to zero significantly before 
the critical temperature is reached, though. 

States with off diagonal long range order by virtue of 77-pairing are defined 
from fermionic operators Cj.|, Cj.j and the fermionic vacuum | 0) as 

L 

Vj — Cj, |Cj- 1 ; ?y := Xl'^J.TCj.i (3.14) 

I*) = 7?^^|0) (3.15) 

These are symmetric states and consequently, their concurrence vanishes in the 
thermodynamic limit due to the monogamy property of pairwise entanglement 
of su{2) degrees of freedom. Consequently, a connection to the order parameter 
of off diagonal long range order 

= (3.16) 



N{L - N) 
L{L-l) " 



>n(l — n) 



(with N, L — > 00 and fixed filling fraction n) can not be established, not even 
for the rescaled concurrence Cr, since 



(see also the analysis for the LMG model in Section [5T^ . Nevertheless, the state 
is entangled, as can be seen from the entropy of entanglement and the geometric 
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measure of entanglement [141] . The latter is tightly connected to the relative 
entropy |140j . Both have been calculated in Ref. [125j and clearly indicate the 
presence of multipartite entanglement. 

3.4.3 Hubbard-type models for interacting fermions 

An interesting class of interacting fermion models is that of Hubbard type mo- 
dels. The Hubbard model [42] is defined by the Hamiltonian 

n = -tY^claCj^a + h.c] + C/^n,,tn,,4 - /iTV (3.17) 
(y) ' i 

where Ci^a, are fermionic operators: {ci.a , c^j = SijSaa'- The coupling 
constant U describes the on-site repulsion, t is the hopping amplitude and fi the 
chemical potential. 

A first study of entanglement in such a system has appeared in [55j for 
the one-dimensional extended Hubbard model for fermions with spin 1/2. The 
extension consists in a nearest neighbor density-density coupling V. Due to the 
conservation of particle number and z-projection of the spin, the local density 
matrix of the system takes the simple form 

P^''' ^z\0){0\+u+\^){U+u-\i){i\+w\tl){n\ (3.18) 
independent of the site number j because of translational symmetry. The broken 




Abbildung 3.7: The local entropy as a function of the on-site coupling U and 
nearest neighbor coupling V. The contour plot below includes the known pha- 
se diagram of the model (full black line). Except the superconducting phase 
transition line, the phase diagram is nicely reproduced. From Ref. [5 5) 

translational invariance in the charge density wave phase has not been taken 
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into account in this work. This does not affect the central resuh but might affect 
the entropy within the charge density wave phase. Except the superconducting 
phase, the phase diagram at half filling (for n — 0) of this model has been 
nicely reproduced by the contour plot of the local entropy (see Fig. 13. 7p , where 
the phase transition lines coincide with its crest. This happens to be an often 
encountered feature of local entropies - also for spin models - as opposed to 
the concurrence for pairwise entanglement whose maxima in general appear at 
a certain distance to quantum critical points and hence are not associated to 
the quantum phase transition. In view of the monogamy of entanglement this 
is evidence for dominant multipartite entanglement in the vicinity of quantum 
phase transitions. 

This analysis clearly points out that the local entropy indicates different 
phase transitions in different ways, essentially depending on whether this quan- 
tity is sensitive to its order parameter or not. Due to the u(l) symmetry of 
the model, the single site reduced density matrix is a functional of occupation 
numbers only. These operators cannot, however, describe order parameters of 
superconductivity or some order parameter of the metal-insulator transition. 
Indeed, the superconducting phase can be predicted if the entropy of entangle- 
ment is calculated for a block of spins, instead of for just a single site [31]. A 
reduced density matrix of at least two sites is necessary for being sensitive to 
superconducting correlations (see also Ref. [69j for a similar result obtained for 
the ionic Hubbard model.) 

Another model studied with this respect is the so called bond-charge exten- 
ded Hubbard model, also known as the Hirsch model, which has originally been 
proposed in the context of high Tc superconductivity [SD] . The Hamiltonian is 

L 

H = U^^ rii^^rii^l - t[l - x{ni^a + '^i+i,-cr)]c|',o-Cj+i,o- + h.c. (3.19) 

■i 

For X = Q the model (|3.19p coincides with the usual Hubbard model p.l7p . In 
phases II and HI (see Fig. 13. 8p there are superconducting correlations due to 77- 
pairing and hence there is multipartite entanglement, as discussed above |125) . 
For X — 1, the model is exactly solvable, and the entanglement of the model has 
been analyzed in Ref. [S] • For general x and n = 1 see Ref. [7J [S] • Besides the 
local entropy of entanglement Si , also the negativity [138] and the quantum mu- 
tual information ^54j have been used and compared for this analysis. While 
Si measures all (pairwise and multipartite) quantum correlations involving this 
specific site, the negativity gives a lower bound for the quantum correlation of 
two specific sites, and the mutual information accounts for pairwise quantum 
and classical correlations. Therefore, this combination of correlation measures 
opens the possibility to decide, what type of correlation is relevant at a quantum 
phase transition. The local entropy is shown in Fig. 13.81 The different phases 
are discriminated by local occupation numbers; consequently, the entropy Si 
bears the information on all the phase diagram except the insulating line IV. 
This is seen from the plot of duSi (with u = U/t) as a function of u and the 
filling fraction n. A comparison of first derivatives respect to y ~ n,u (depen- 
ding on the phase transition) of all three correlation measures reveals common 
singularities for dySi and dyTij only for the transitions II-III and II-IV; fur- 
thermore, it is found that the range R of the concurrence diverges[H| at both 
transitions. These facts allow to judge the transitions II-III and II-IV (at n = 1 
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Abbildung 3.8: Left panel: The ground state phase diagram of the Hirsch model 
at a; = 1. Empty, slashed and full circles indicate the presence of empty, singly 
and doubly occupied sites, respectively. Right panel: Except the insulating line 
IV, the phase diagram is nicely reproduced by duSi. [From [5]] 

and arbitrary x) as governed by pairwise entanglement. For the transitions II-I 
and II-I' instead, multipartite entanglement is relevant, accompanied by a finite 
range of the concurrence. A similar behavior was encountered for non-critical 
spin models where the divergence of R was accompanied by the emergence of a 
fully factorized ground state. Here, R diverges close to QPT. 

In order to detect the charge transition II-IV a,t n — 1 and any x, dxSi has 
been calculated by means of DMRG . Its singularities allowed to accurately 
determine the charge gap as a function of the bond-charge coupling x. 

We finish this selection of results for itinerant fermion systems with the Hub- 
bard model in a magnetic field. Also here, the local entropy Si has been looked 
at in order to analyze its entanglement. As in the examples before. Si indicates 
the second order phase transitions in terms of divergences of its derivatives dhSi 
and d^Si, respectively. Indeed, it has been demonstrated that duSi and df^Si can 
be expressed in terms of spin and charge susceptibilities [67j . This bridges ex- 
plicitly between the standard method in condensed matter physics for studying 
phase transitions and the approach from quantum information theory. 

Summarizing, the body of work developed so far suggests the conclusion 
that local entropies can detect QPTs in systems of itinerant fermions, particu- 
larly if the transition itself is well predicted by a mean field approach for local 
observables of the model. In the described cases, translational invariance leads 
to predictions independent of the site, the local entropy is calculated for. In 
absence of this symmetry, it might prove useful to average over the sites; the 
resulting measure is then equivalent to the Q-measure [S2] . 

Though there certainly are transitions with dominant features in the pairwise 
entanglement, also here the generic case indicates the dominance of multipartite 
quantum correlations. 

3.4.4 A remark on entanglement of particles 

There is little work which uses measures for indistinguishable particle entan- 
glement (see Section l2.5p . particularly regarding the use of the fermionic con- 
currence, giving account for the possibility of double occupancy (with internal 
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degree of freedom). The entanglement of particles (see Section [^33]) and its dif- 
ference from the usual spin entanglement has been worked out in [37j ; starting 
with very small systems as two spinless fermions on four lattice sites and the 
Hubbard dimer, and then for the tight binding model in one spatial dimension, 
the results are compare with previous results for the spin entanglement in |124j . 

For the Hubbard dimer (a two-site Hubbard model), the authors compare 
with the results for the entanglement measured by the local von Neumann entro- 
py without superselection rule for the local particle numbers |148j . Whereas the 
latter signals decreasing entanglement in the ground state with increasing U /t, 
the entanglement of particles increases [37]. This demonstrates that imposing 
superselection rules may lead to qualitatively different behavior of the entangle- 
ment. Interestingly, an increase with U /t is observed also for the entanglement 
of modes (without imposing superselection rules jl. 

We would like to finish this section with the notice of a recent proposal of 
an experiment in order to decide whether 'entanglement' merely due to the sta- 
tistics of the indistinguishable particles can be useful for quantum information 
processing [27] . 



^For the extended Hubbard dimer, which is defined as the two site extended Hubbard 
model, see I35| 
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Kapitel 4 

Multipartite entanglement: 
quantification and 
classification 



One conclusion from the preceding chapter is that multipartite entanglement 
plays an important role in condensed matter systems, in particular close to a 
quantum phase transition. In this chapter we will present a whole variety of 
quantities that have been proposed for the scope of quantifying and eventually 
classifying multipartite entanglement. 

4.1 A zoo of multipartite entanglement measu- 
res 

Both the classification of entanglement and its quantification are at a prelimina- 
ry stage even for distinguishable particles (see however |3S1 [SH 11301 [HH 
inn HD] and references therein). This uncertainty is responsible for the vast 
amount of suggested quantities for the analysis of multipartite entanglement 
in many-particle systems and in particularly interesting wave functions. It has 
already been mentioned that several quantities applied in the previous sections 
are useful as indicators for multipartite entanglement when the whole system is 
in a pure state; then the cumbersome convex-roof construction is not needed. 
For non-degenerate ground states of model Hamiltonians, this requirement is 
met at zero temperature T = 0. The entropy of entanglement is an example 
for such a quantity and several works use multipartite measures constructed 
from and related to it (see e.g. [501 [5^ [T^ 11151 [551 [75] ). These measures are of 
'collective' nature - in contrast to 'selective' measures - in the sense that they 
give indication on a global correlation without discerning among the different 
entanglement classes encoded in the state of the system. The advantage of these 
measures is that they are easily computed. Their disadvantage is founded in an 
ambiguity concerning their choice: as soon as at least two different entanglement 
classes are measured, an infinite variety of inequivalent measures does exist. As 
an example, infinitely many proposals could be generated from [1151 1551 15H [75] 
by substituting the standard entropy with another mixedness measure (e.g. the 
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one-parameter family of Renyi entropies, etc.). Since these new measures will 
induce a different ordering in the space of entangled states, some results will 
deviate. The main problem would then consist in extracting the information 
about different entanglement classes from a vast collection of results. In this 
context, the authors of Ref. [43l [44] suggest the analysis of a distribution of 
purities for different bipartitions. If such a distribution is sufficiently regular, 
its average and variance might prove characteristic for the global entanglement 
in the system. The application of this analysis to the one-dimensional quantum 
Ising model in a transverse field revealed a well behaved distribution function, 
whose average and second moment are good indicators of the quantum phase 
transition |31j : at the quantum critical point both the average and the standard 
deviation exhibit a peak that becomes more pronounced as the number of qubits 
is increased. 

Another pragmatic way of quantifying entanglement in a collective way is 
represented by the geometric measure of entanglement |141j . It quantifies the 
entanglement of a (pure) state through the minimal distance of the state from 
the set of (pure) product states [1271 1141j 

i?<,(vE') = -log2max|(*|$) P (4.1) 

where the maximum is on all product states $. As discussed in detail in |141| . 
the previous definition is an entanglement monotone. It is zero for separable 
states and rises up to unity for e.g. the maximally entangled n-particle GHZ 
states. The difficult task in its evaluation is the maximization over all possible 
separable states and of course the convex roof extension to mixed states. 

A different approach was pursued in [57] (see also |116] ) where different 
bounds on the average energy of a given system are obtained for different ty- 
pes of n-particle quantum correlated states. A violation of these bounds then 
implies the presence of multipartite entanglement in the system. The starting 
point of Giihne et al. are the n- separability and k-producibility which admit to 
discriminate particular types of n-particle correlations present in the system. A 
pure state | '0) of a quantum systems of N parties is said to be n-separable if it 
is possible to find a partition of the system for which | ip) = \4>i)\4>2) ■ ■ ■ {(f'n)- A 
pure state | ip) can be produced by k-party entanglement (i.e. it is k-producible) 
if we can write \ ip) = \(j)i)\(j)2) ■ ■ ■ \(j)m) where the \(f>i) are states of maximally k 
parties; by definition m > N/k. It implies that it is sufficient to generate spe- 
cific k-party entanglement to construct the desired state. Both these indicators 
for multipartite entanglement are collective, since they are based on the facto- 
rizability of a given many particle state into smaller parts, k-separability and 
-producibility both can not discriminate the different k-particle entanglement 
classes (as e.g. the k-particle W-states and different k-particle graph states [58j . 
like the GHZ state). 

Another approach pursued is the generalization of the celebrated concur- 
rence. For the quantification of pairwise entanglement in higher dimensional 
local Hilbert spaces, the concept of concurrence vectors has been formula- 
ted [ini [H] besides the I-concurrence |108j ; the length of the concurrence vector 
has proved equivalent to the I-concurrence |146j . Also for multipartite systems 
of qubits the concurrence vector concept has been proposed 2 . In the multi- 
partite setting however this means to apply the pure state concurrence formula 
to a mixed two-site reduced density matrix. It will coincide with the true con- 
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currence if and only if the eigenbasis of the density matrices accidentally are an 
optimal decomposition. Therefore, the concurrence vector in this case is not a 
vector whose entries are the concurrences, and it is at least not obvious whether 
this proposal is an entanglement monotone. 

The n-tangle is a straightforward extension of the concurrence to multipar- 
tite states as the overlap of the state with its time- reversed |144| . It vanishes 
identically for an odd number of qubits, but an entanglement monotone is ob- 
tained for an even number of qubits. Due to its factorizing structure, it detects 
products of even-site entangled states in addition to certain genuine multipar- 
tite entangled states: it detects the multipartite GHZ or cat state, but not for 
example the four qubit cluster state. Therefore, also the n-tangle is a collective 
measure. 

4.2 Measures for genuine multipartite entangle- 
ment 

The counterpart to collective entanglement measures are selective measures for 
each different multipartite entanglement class separately. A first multipartite 
example beyond pairwise qubit entanglement has been derived from the concur- 
rence and the (one-)tangle: the 3-tangle [30] . It is a measure for genuine tripartite 
entanglement that discriminates also from pairwise entanglement distributed all 
over the chain, as is the case for the W state | VF) ~ 1 100) + \ 010) + \ 001). The 
3-tangle coincides with the 3-dimensional hyperdeterminant for two-dimensional 
local vector spaces [H3] , i.e. a generalized determinant form for 3x2 matrices over 
C. It originated the insight that SL{2, C) invariance rather than SU (2) invarian- 
ce leads to a classification of genuine multipartite entanglement: the 3-tangle is 
indeed the only SL{2, C)®^ invariant, and only a single entanglement class with 
respect to SLOCC does exists. The task of finding SLOCC-class selective 
measures of qubit entanglement hence is reduced to finding local SL(2,<C) in- 
variant operators, from which global 5^(2, C)®" invariants can be constructed. 
The invariance group should be extended to including the symmetric group 5*^ 
for g-partite systems [93|3. The genuine multipartite entangled states belong to 
the non-zero SLOCC dasB 

There are standard formalisms in the well developed field of invariant theory 
for the construction of all SL{2, C)**" invariants (see Ref. |22j and references 
therein, but also Ref. [129j for the contraction method with the invariant tensor, 
which for qubits or spins-1/2 is i(Jy). However, the construction of a complete set 
is very cumbersome already for five qubits |79j and the so constructed invariants 
will typically be non-zero also for certain factorizing states, and consequently 
will they be collective measures. Furthermore would we desire to walk in the 
steps of Wootters' concurrence and to write the invariants as expectation values 
of some tangle operator. The minimal requirement would then be that the mea- 
sure of a somehow factorizing state should be zero. Interestingly, it seems that 

^This amounts to including the permutation of qubits to the realm of "local" operations; 
this bigger symmetry group has been termed SLi^c* in Ref. 29 . 

^ The complementary zero SLOCC class is made of those states that vanish after - maybe 
infinitely many - suitable SL{2, C) operations; the non-zero SLOCC class is the complement 
of the zero SLOCC class. 
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this 'minimal' requirement already implies the necessary invariance, at least for 
qubits [53 [Ml- 

The special case of factorization into a single qubit and the remaining (n-1) 
qubits leads to the central quality of the local operators we are searching for: 
its expectation value should vanish for all single qubit states | V') 

(V'|O|V) = (4.2) 

Such an operator can impossibly be linear, and antilinear operators have to be 
studied [53]. The unique antihnear operator satisfying Eq. (|4.2p is ay€, as used 
for the definition of the pure state concurrence, where £ is the complex conjuga- 
tion in the eigenbasis of <7z ■ It is suggested from the described shortcomings of 
the n-tangle as a measure for genuine multipartite entanglement measure jl44] 
that at least one further operator besides ay will be needed. Such an additional 
operator is found when the requirement (j4.2p is extended to the two-fold copy 
of single qubit states 

(V I . (V- 1 O I V-) • I V-) = (4.3) 

The requirement to being trace-norm orthogonal to ayuay€ leads to the unique 
(up to normalization) operator 

3 

^1=0 

where 5^ : {go, gi, 92, 93) — IjO, 1) and cq = l2- Both antilinear operators 
are invariant under 5*^(2, C) operations on the qubit. 

One possible way of expressing the 3-tangle in terms of the above mentioned 
pair of antilinear operators is 

= \ (-0* I • (V'* «) CT,. CTA «) cr^ «i cr'' (g) cr^ I 7/)) • I • (4.5) 

This form has the minimal possible multi-linearity and it is evidently per- 
mutation invariant. Unfortunately, the convex roof construction presented in 
Refs. |145l I123j for bilinear quantities can not in general be extended to higher 
multi-linearity. This is very much related to the absence of a tilde- orthogonal 
basis |145| that exists for bi- antilinear constructions (see |123| ). 

A first success in direction of an understanding of convex roofs for multipar- 
tite entanglement measures is the analytic convex roof construction found for 
rank two- mixtures of translational invariant GHZ and orthogonal W states [75] . 
The result is best expressed in the Bloch sphere of both orthogonal states (left 
panel of Fig. 14.11) . The right panel of Fig. 14.11 shows the convex roofs of the 
three entanglement measures tangle(solid line), concurrence (dashed line) and 
3-tangle (dotted line), which enter the monogamy relation of Coffman et at. 
The curious result is that an interval with vanishing concurrence and 3-tangle 
appears, whereas the tangle is always positive. On the background that only 
two classes of entanglement exist for three qubits [38] , this result might appear 
contradictory at first sight. That this is not the case, is due to a subtlety of con- 
vex roof measures, when different entanglement classes are considered (here the 
concurrence - pairwise - and the 3-tangle). The explanation is that the optimal 
decompositions leading to the respective convex roofs, will in general not be 
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GHZ (p=l) 




Abbildung 4.1: Left panel: Bloch sphere for the two-dimensional space spanned 
by the GHZ state and the W state. The simplex SO contains all states with zero 
three-tangle. The leaves between pO and pi represent sets of constant three- 
tangle, and in the simplex SI the three-tangle is afhne. The corners of the 
simplices constitute the optimal decomposition. Right panel: Plot of the convex 
roofs of the tangle (solid line), concurrence (dashed line) and 3-tangle (dotted 
line). There is an interval with entangled states without concurrence or three- 
tangle. [From [75] ] 



compatible, in the sense that both will not emerge from a single optimal decom- 
position of the whole density matrix. The physical interpretation of this curious 
finding is very simple: in the interval of vanishing concurrence and 3-tangle, the 
density matrix can be built from pure states that have either concurrence or 
3-tangle (or both). 

With the two independent local S'L(2, C)-invariant operators (combs) (|4.2p 
and (j4.3p at hand, entanglement measures for genuine multipartite entanglement 
can now be constructed for an arbitrary number of qubits. The astonishing 
news from the four qubit case is that three inequivalent measures for SLOCC 
entanglement classes do exist 

(af^a^ffyay) • (a'^ayaxay) • {(Jya^a^ay) (4.6) 
&{a^a,^(7yay) • {a^(Jy(J\ay) • 

• {oycy^ (TyGr) • {ayayG^a'^) (4-7) 

• (aPaya'^ay) • {(JyCFxcr^ay) • (aya^a'^ay) (4.8) 

where 6 indicates the symmetrization with respect to the permutation group 
S'4 on four qubits. This leads to seven inequivalent representatives of SLOCC 
entanglement classes for four qubits H These include the GHZ state (|4.9p . the 
celebrated cluster state ()4.10p (three states connected by means of ^4 permuta- 
tions) and a third state (|4.11|) . which in contrast to the former two has non-zero 

^The seven representatives are characterized by their different characteristics on the three 
measures (see Table : three each with one or two measures being non-zero plus one where 
all three measures are non-zero. 
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4-qubit hyperdeterminant. 

|$2> = -^(10000) + I nil)) (4.9) 
I $4) = i(| nil) + 1 1100) + 10010) + 10001)) (4.10) 

I $5) ^(\/2 1 nil) + 11000) + 10100) + 10010) + 10001)) (4.11) 

Furthermore, representatives of g — 1 inequivalent non-zero SLOCC classes [551 
I129j are known for q qubits [Mj. Explicit measures for up to six qubits and a 
prescription how general g-qubit entanglement measures are constructed have 
been given in Ref. [H] ; their characteristics on g — 1 maximally entangled states 
demonstrates that the states belong to different non-zero SLOCC classes (see 
table 1121). 

Tabelle 4.1: Filter values for maximally entangled states; the length is the num- 
ber of Fock-elements in their canonical form. The 5- and 6-qubit entanglement 
measures are explicitly given in Ref. 94J together with the maximally entangled 
states they have been evaluated on here. 
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Summarizing, a method has been constructed that permits to construct 
measures for genuine multipartite entanglement, i.e. for the non-zero SLOCC 
class of entangled states. Eventually necessary symmetrization with respect to 
the corresponding symmetric group leads to generators of the ideal of SL*- 
invariants that vanish on arbitrary product states. The generators of this 
ideal can also be constructed from the invariants known from invariant theo- 
ry [711 [75] . It would be highly desirable to bridge between both approaches. 
For four qubits, this bridge has already been established by the author in colla- 
boration with D. Z. Dokovicto the extent that the three 4-qubit entanglement 
measures have been expressed in terms of three generators of that ideal but con- 
structed from the standard invariants from [751 US]. The results will be presented 
in a forthcoming publicatiorQ. 

*Note added after the habilitation procedure: this work is now available on the ar- 
Xiv:0804.1661 
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4.3 Experimental access to genuine multipartite 
ent anglement 

The invariants under local SL(2,<C) operations are most naturally expressed 
in terms of expectation values of antilinear operators. However, in the labora- 
tory physical observables are measured, which are linear operators. So, unless 
the experiment is performed with quantum states with real coefficients only, a 
translation of the above multipartite entanglement measures into spin correla- 
tion functions will be necessary. Such a one-to-one translation indeed exists. In 
order to see how the modulus of antilinear expectation values can be transported 
into linear expectation values, we write 



{i;\{d<E> ... )\r)---{r\{d<E>...)\i^)... 

= (^i.(^*i....[(o.onpMO. 



(4.12) 
D...|V)-|^*>---- 



where P^^i' is the permutation operator acting on the corresponding different 
copies of the same qubilH and the product extends over all qubits and various 
copies such that exclusively linear expectation values occur. As a consequence, 
the linear operator corresponding to an antilinear operator O is then 



£[0] := (l.e:)(O.C')P 



(4.13) 



It is clear that the transformation in the opposite direction works the same way. 
In particular we find 



diag {1,-1, -1,-1} 



(4.14) 
(4.15) 



It is worth noticing that Af^i, is precisely the Minkowski metric. 



(4.16) 
(4.17) 



-2 -2 \ 





\ -2 -2 y 



, J = 



/ 2 2 \ 

110 

110 

V 2 2 / 



(4.18) 



This is a straight forward prescription for expressing the genuine multipartite 
entanglement measures presented in Ref. [93, 94J in terms of expectation values 
of linear operators. The details of this transcription including an analysis of 
weakly mixed states will be presented in a forthcoming publication. 



= • (j^, using Einstein sum convention 
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Kapitel 5 

Conclusions & Outlook 



In conclusion, quantum entanglement is a fascinating phenomenon that natural- 
ly occurs in the eigenstates of interacting Hamiltonians. It is certainly the inter- 
pretation of this entanglement as a resource for quantum information processing 
that nowadays is mainly responsible for its growing relevance and consideration 
in e.g. condensed matter physics. Nevertheless did measures for certain classes 
of entanglement also prove to be useful tools for the indication of quantum pha- 
se transitions if the corresponding measure (e.g. some entanglement entropy) 
incorporates the relevant order parameter. Though there is only a single class- 
selective entanglement measure we can compute for mixed states, namely the 
concurrence for pairwise entanglement of qubits, the Coffman-Kimdu-Wootters 
relation makes the local von Neumann entropy to an implicit indicator of multi- 
partite entanglement. The analysis of ground state entanglement of many model 
Hamiltonians gives strong support to the conclusion that multipartite entangle- 
ment becomes dominant in quantum Hamiltonians close to their critical point. 
But it is also rather multipartite than only pairwise entanglement, that is crea- 
ted by Hamiltonian evolution from disentangled initial states. A deeper analysis 
is needed here in particular in the direction of more conclusive measures rather 
than analyzing more and more different Hamiltonians with always the same ca- 
nonical set of measures. On the other hand, since the presence of multipartite 
entanglement seems to be the generic case rather than the exception, it might 
be the weakly or even disentangled states, as for the factorizing field in spin 
chains, which are the precursors of drastic changes and complex phenomena. 
In any case must the understanding of entanglement and its subdivision into 
distinct SLOCC classes be improved substantially before a possible connecti- 
on between entanglement and quantum phase transitions can be established. A 
step into this direction is to realize that the relevant symmetry group for entan- 
glement classification is the local 5L(2,C), including the permutation group of 
the qubits, and that entanglement monotones for this invariancc groiip can be 
built from local invariant and antilinear operators with zero expectation values 
on all the local Hilbert space - combs. Two such operators have been found 
recently for qubits. They permit the construction of a complete set of measures 
for genuine multipartite entanglement provably up to 4 qubit systems. The main 
achievement lies in that the combs admit a straight forward construction of the 
ideal of SL{2, C)'^^-invariants whose elements vanish on all product states and 
for arbitrary number of qubits N; this is a problem which already for five qubits 
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creates severe difficulties for standard metliods from invariant theory, where it 
seems hopeless to handle them for systems larger than this. A one-to-one corre- 
spondence of the antilinear invariants to correlation function finally even gives 
a prescription for the laboratory in order to extract these quantities from spin 
correlation functions. Yet, this is rather the very beginning of research in this 
direction, since there is little understanding of even very elementary questions. 
Although a set of maximally entangled states is given for an arbitrary num- 
ber of qubits, the ambiguity of such a set is an issue: it is not clear from the 
beginning, whether a complete set of maximally entangled states can possib- 
ly be constructed from them as elements of kind of a basis for such states. A 
characterization of maximally entangled states might be thinkable in terms of 
the analogue of a tilde- orthogonal basis (see [145] ). Completeness of the ideal 
is another big issue which would have an important impact also in the field of 
invariant theory. A further major challenge consists in the extension of the pure 
state measures to mixed states, which at the end is necessary in order to express 
the entanglement measure in terms of correlation functions of a practically infi- 
nitely large condensed matter system. Also here an important step forward has 
been done on the mixture of three qubit GHZ and W states, but the ultimate 
goal would be an analytic procedure for the convex roof extension along the 
lines of the concurrence. Already now the accumulated knowledge about opti- 
mal decompositions is of substantial help for numerical procedures and leads to 
nontrivial lower bounds for multipartite entanglement in mixed states, but the 
major part of the puzzle pieces is yet to be found. A consequent extension of 
the method of combs would consist in their identification for local Hilbert space 
dimension larger than two. Here the problem is more complex already at the 
starting point, because there are no bi-antilinear SL{2S + 1, C) invariant combs 
for spin S > Last but not least we mention that once, genuine multipartite 
entanglement measures are known, they should be used for testing and ana- 
lyzing alternative approaches to the detection of entanglement. Witnesses are 
one relevant example of a complementary approach, which is motivated from 
the question for separability of a mixed quantum state rather than from the 
classification of entanglement. It might be just their complementarity bearing 
the key towards many open questions in the theory of entanglement. 
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